Continuous-wave spatial quantum correlations of light induced by multiple scattering 



(N 

o 

(N 



Stephan Smolka^'B Johan R. Ott, 1 Alexander Huck, 2 Ulrik L. Andersen, 2 and Peter Lodahl 3 '0 

1 DTU Fotomk, Department of Photonics Engineering, 
Technical University of Denmark, Building 345V, 2800 Kgs. Lyngby, Denmark 
2 DTU Physics, Department of Physics, Technical University of Denmark, Building 309, 2800 Kgs. Lyngby, Denmark 



Niels Bohr Institute, University of Copenhagen, 



amsvej 17, Dk-2100 Copenhagen, Denmark 



We present theoretical and experimental results on spatial quantum correlations induced by mul- 
tiple scattering of nonclassical light. A continuous mode quantum theory is derived that enables 
determining the spatial quantum correlation function from the fluctuations of the total transmit- 
tance and reflectance. Utilizing frequency-resolved quantum noise measurements, we observe that 
the strength of the spatial quantum correlation function can be controlled by changing the quantum 
state of an incident bright squeezed-light source. Our results are found to be in excellent agreement 
with the developed theory and form a basis for future research on, e.g., quantum interference of 
multiple quantum states in a multiple scattering medium. 
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I. INTRODUCTION 

The light transport through an inhomogeneous ran- 
dom medium is determined by a multiple scattering pro- 
cess whereby interference between different light paths 
induces a complex intensity speckle pattern. Mesoscopic 
intensity fluctuations that survive averaging over all con- 
figurations of disorder give rise to phenomena such as 
enhanced backscattering and Anderson localization of 
light* 1 - In order to characterize such a disordered medium, 
it is essential to investigate spectral or temporal inten- 
sity correlations in the multiply scattered light, 2-4 In the 
diffusive regime, the light intensities of different spatial 
directions (i.e. independent speckles) are uncorrelated. 
In this case interference effects disappear after ensemble 
averaging and the transport is well described by diffusion 
theory. 

While the classical aspects of multiple scattering have 
been studied extensively, the quantum nature of light in 
a multiple scattering setting has been addressed only re- 
cently. Triggered by theoretical studies on quantum fluc- 
tuations of light in absorbing and amplifying multiple 
scattering media^£ recent experiments have addressed 
the quantum noise of multiple scattered lighter— and spa- 
tial photon correlations J£r— Furthermore, first studies of 
photon emission and cavity quantum electrodynamics in 
disordered media have been carried out 1 ^— . It was re- 
cently predicted that quantum interference and entangle- 
ment can be induced by combining two or more quantum 
states in a multiple scattering mediuni ) 18 i 19 which could 
be of potential interest for applications in quantum in- 
formation processing. In the diffusive regime of multi- 
ple scattering classical intensity correlations vanish and 
many quantum optical phenomena do not survive ensem- 
ble averaging including polarization entanglement^ and 
quadrature squeezing^ In contrast, it has been proposed 
that nonclassical photon fluctuations can be detected in 
the multiple scattered light even after ensemble averag- 
ing, which results in a new type of correlations that is 
of quantum origin. 10 These spatial quantum correlations 



have been reported recently based on total transmission 
quantum noise measurements, thereby entering for the 
first time experimentally the genuine quantum optical 
regime of multiple scattering^ 

In this report, we present a detailed theoretical and 
experimental analysis of spatial quantum correlations in 
transmission and reflection measurements. We develop 
in Sec. II a continuous mode quantum theory in order 
to relate the spatial quantum correlation function to the 
photon fluctuations of the light source. By exploiting 
quantum noise measurements, the variance of the pho- 
ton fluctuations of the multiple scattered light is studied 
in Sec. Ill from which the strength of the spatial quan- 
tum correlation function is determined. Finally, it is ex- 
perimentally shown how the spatial quantum correlation 
is tuned by controlling the quantum state of the light 



II. THEORY 



The quantum nature of light is revealed, e.g., by in- 
vestigating photon fluctuations and spatial, temporal, or 
spectral quantum correlations among photons. In the 
experiment of concern here, the quantum properties are 
identified through the photon fluctuations for which sub- 
Poissonian fluctuations is a sign of non-classicality^ A 
schematic illustration of the experimental setup is dis- 
played in Fig. [T] a bright squeezed state of light is sent 
through a multiple scattering medium, the total trans- 
mitted or reflected light is collected, and the associated 
photon fluctuations are recorded. In the following, we 
introduce a continuous mode quantum theory that is di- 
rectly applicable to the quantum noise measurements dis- 
cussed in the present work. 

We define the spectral density function at frequency 
oj as the Fourier transformation of the autocorrelation 
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FIG. 1. Illustration of multiple scattering leading to spatial 
quantum correlations. The nonclassical light source that is 
used in the experiment is created by overlapping a squeezed 
vacuum light beam (SQZ) and a coherent light beam (C) on 
a beam splitter with two input ports, 1 and 2, respectively. 
The relative phase, A(f>, between C and SQZ can be tuned 
continuously. The resultant light beam in direction a is in- 
cident onto a medium with randomly distributed scatterers. 
Light is split into a multitude of different trajectories and the 
number of photons exiting the medium in a specific direction, 
b, can be correlated with the number of photons in another 
direction, b' , to a degree dependent on the quantum state 
of the illuminating light source. In the diffusive regime the 
spatial quantum correlation function can be determined from 
total transmission quantum noise measurements. A similar 
experimental scheme is applied for reflection measurements. 



function^ 



where 



r(r) = (n(t)n(t + r)> - (h(t))(h(t + r)>, (2) 



with n(t) being the operator describing the photon flux. 
The variance of the photon fluctuations is related to 
the noise spectrum through the inverse of Eq. (JXJ) , i.e. 
Ah 2 (t) = J^^dujS^). The spectral density is studied 
within a frequency region uj± = co s ± Slu/2 
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that is proportional to the experimentally measured 
noise power. 

The randomly disordered medium is treated as N spa- 
tially distinct optical input modes that are intrisically 
connected to N output modes via multiple scattering of 
light. In order to relate the quantum properties of the 
multiple scattered light to the quantum state of the inci- 
dent light source, the spatial output mode b is related to 
all input modes a' through 



a b (t) 



, S a 'b 



a a > (t), 



(4) 



where a b (t) denotes the time-dependent annihilation op- 
erator and s a 'b is the complex electric field amplitude 
scattering coefficient corresponding to either a reflec- 
tion (r ab ) or transmission (t a b) channel. The corre- 
sponding intensity transmission (reflection) coefficient is 
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The annihilation and creation 



operators obey the commutation relation [at(t), a b , (t')] = 
Sb.b' 8(t — t'), with 5(t — t') and 8 b ^ being the Dirac 
and Kronecker delta functions, respectively. The av- 
erage photon flux is defined as (hb{t)) = (al(t)a,b(t)) 
while the variance in the photon fluctuations is An 2 (t) = 
(n 2 (t)) — (nfe(i)) 2 , where (• • • ) describes the quantum 
mechanical expectation value. In the following we fo- 
cus on the setting where the light beam is incident onto 
the medium through a single direction, a, thus the aver- 
age photon flux for all other input modes equals zero. 
The mean photon flux of the total transmission can 
be obtained by summing over all transmitted output 
modes, i.e. (hr(t)) = J2b Tab(a a (t)a a (t)) , where the in- 
tensity transmission coefficients are random variables de- 
termined by the individual realization of disorder. Pre- 
dictable physical variables are obtained after ensemble- 
averaging over all configurations of disorder, e.g., the to- 
tal transmitted photon flux is given by 



(nr(*)> =T<n(t)>, 



(5) 



where T = J2 b Tab is the total intensity transmission for 
the flux of photons (n(t)} incident in mode a, and the 
bar denotes the ensemble average over all configurations 
of disorder.— A similar result can be found for the total 
reflection, i.e. (n,R(t)) = R(n(t)) where R = J2b^- ab - 
In order to determine the fluctuations of the total trans- 
mission we need to evaluate the second moment of the 
photon number operator 



(n T (t)n T (t + t)) = {h b {t)h b ,{t + T)) 



b,b' 

=T*(h(t)n(t + t)) + {T- W)(n(t))6(T). 

(6) 

The non- vanishing contributions to the sum in Eq. for 
different output modes b and b' stem from the fact that 
different spatial parts of the volume speckle patterns are 
quantum correlated. In an experiment the total collec- 
tion and detection efficiency is non-unity, which is ac- 
counted for by substituting T — > t\tT . The photon fluc- 
tuations of the total transmission, (An^(w s ,&)), can 
be related to the spectral density function of the inci- 
dent light source integrated over the measured frequency 
range from Eqs. <jT][3]) 
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where we have used T 2 = T that holds in the diffusive 
regime. 4 Evaluating the last expression for a coherent 
state, yields (An 2 -, T (u a , Soj)) — rfT(hc(t)) 5uj/2it. 
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We introduce the Fano factor 

F{u,6u) = 2 — j-^r. 

(An£(w,<>w)) 



(8) 



that gauges the ratio between the variance in the 
photon fluctuations of an unknown quantum state, 
(An 2 (oj s , Sui)), and of a coherent state, (Anp(w s , 
having the same mean photon flux. Importantly, the 
Fano factor can be directly measured without determin- 
ing the proportionality factor between the noise power 
recorded with the electrical spectrum analyzer and the 
photon number variance. In the experiment the optical 
state is created by overlapping a squeezed- vacuum state 
with a bright-coherent state, cf. Fig. [TJ Since the re- 
sultant bright-squeezed state has a large coherent ampli- 
tude, the average photon flux is to a very good approxi- 
mation equal to the average photon flux of the coherent 
state and Eq. ([7]) can be rewritten as 

^T 2 ((An 2 K,fe))-("W)^/2vr) , 1 
F t (uJs,ouj) = ==ttt-t, : — - V 1 



r) T T(h(t))6uj/2ir 
--l-r lT T[l- F{cj s ,Scj)}. 



(9) 



A similar expression holds for the total reflection, i.e., 
Fr(lj s ,6ui) = 1 - tjrR [1 - F(u s ,5u)] with r] R being the 
corresponding detection efficiency We emphasize the im- 
portance of the last expressions. In Ref. [21| it was pre- 
dicted that any nonclassical features in the electric field 
quadrature amplitudes vanish in the multiple scattering 
process after ensemble averaging. In contrast, Eq. © 
shows that for a quantum state with nonclassical photon 
fluctuations also the ensemble-averaged transmitted and 
reflected photon fluctuations exhibit nonclassical photon 
fluctuations after multiple scattering. 

We now introduce the spatial quantum correlation 
function, which is the quantity extracted in the exper- 
iment. It is defined as 

C b Q b ,(w s ,<M = — / dr I dujC bb ,(T)e^ T , (10) 
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for b ^b' where 



= (n b (t)n v (t + r)) (u) 
(n b (t)) x {h b ,{t + T)) 

The correlation function gauges the photon correlations 
between two different output modes b and b' . Starting 
from Eq. (j4]) we derive 



[F(u s ,8lj) - 1] Sw-^- 
(n(t)> 2tt fcb " 



(12) 



by using that (n(t)) w (n(t))c — 2tt/Soj (Aiip(w s , Sw)), 
Here C bb , = T ab T ab > / (T ab x T ab > ) is the classical speckle 



correlation function that is induced by multiple scatter- 

9-4 

mgi— . 

From Eq. (fT2")l it can be seen that different optical 
modes b and b' are positively correlated, i.e., the photons 

show spatial anti- bunching (CR, (lj s , Su)) > 0) when the 
light source exhibits classical photon fluctuations. For a 
light source that exhibits nonclassical photon fluctuations 
(F(u> s ,8u>) < 1), we enter the purely quantum optical 
regime where photons propagating two different direc- 
tions have negative spatial correlations (C® b ,(ui s ,5uj) < 
0). The spatial quantum correlation function can be ex- 
pressed in terms of the ensemble averaged total trans- 
mitted (reflected) Fano factor that is directly measured 
in the experiment. This is seen by substituting Eq. (|9]) 
that holds in the diffusive regime {C bb , — 1) into Eq. (fT2l) . 
leading to the expression 



Su) 



1] 5uj 
2^ 



tE- (13) 



Thus, the spatial quantum correlation function can be ex- 
tracted by recording in addition to the Fano factor also 
the photon flux and the detection efficiency for the to- 
tal transmission or reflection from the multiple scatter- 
ing medium, as schematically shown in Fig. [TJ The mean 
photon flux can be measured by recording the light power 
of the light source, P, divided by the energy of a photon 
at the optical frequency, loq, such that (n(t)) = P/Hujq, 
where h is the Planck constant divided by 2 n. Our 
method provides an experimentally simple way of ex- 
tracting two-channel speckle correlations in comparison 
with a direct measurement based on correlation measure- 
ments with two single photon counting detectors. 13 



III. EXPERIMENT 

In the following the performed experiment is discussed. 
First, the experimental setup and the characteristics of 
the used light source are described in Sec. III. A. A de- 
tailed characterization of the multiple scattering sample 
is carried out in Sec. III.B. In Sec. III.C we present the 
experimental results on the measured multiple scattered 
photon fluctuations and show how the magnitude of the 
spatial quantum correlations can be modified. 



A. The squeezed light source 

The experimental setup to study the photon fluctua- 
tions of multiply scattered light is schematically shown 
in Fig. [2ja,. The squeezed vacuum state is generated with 
a bow-tie shaped optical parametric oscillator operating 
below threshold . 24 ' 25 As a nonlinear medium, a type III 
periodically poled nonlinear crystal was placed inside the 
optical parametric oscillator and pumped with a contin- 
uous wave laser (Diabolo, Innolight) at a wavelength of 
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FIG. 2. (color online), a, Sketch of the experimental setup. 
MC: mode cleaner, PPKTP: periodically poled potassium ti- 
tanyl phosphate nonlinear crystal as part of an optical para- 
metric oscillator (gray shaded box), SQZ: vacuum squeezed 
light beam, C: coherent light beam with a variable optical 
phase, Atp, relative to the squeezed vacuum, BS: beam split- 
ter, FM: flip mirror, S: multiple scattering sample, SA: spec- 
trum analyzer, D: detector. The detection detector, D2, has 
a bandwidth of 315 kHz. The multiple scattered light is col- 
lected with a microscope objective featuring a numerical aper- 
ture of NA=0.63. In reflection geometry the incident light 
illuminates the sample under an angle a — 69°. A detailed 
explanation of the setup can be found in the main text, b, 
Measured Fano factor, F(co s , Sco), at an elecrical sideband fre- 
quency of oj s = 3.93 MHz of the bright squeezed light source 
without insertion of the multiple scattering sample (black 
circles) depending on the phase A<f> of the coherent beam, 
as schematically shown in the inset. The measured classi- 
cal limit is plotted with gray triangles and is obtained by 
blocking the squeezed beam. As the resolution bandwidth we 
used 8to/2n — 300 kHz. The uncertainty in the Fano factor is 
estimated over 18 full periods of A(j> — . . . 367r. c, Recon- 
structed Wigner function of the squeezed vacuum state (SQZ) 
produced by the optical parametric oscillator. The Wigner 
function is obtained from time resolved measurements using 
a homodyne setup with a 50:50 beam splitter (see inset) and 
quantum state reconstruction. 



Xp — 532 nm. A second-order nonlinear process generates 
photon pairs that are in resonance with the cavity and are 
centered around Ao = 1064 nm. A fraction of the cavity 
mode is coupled out by a partly transmitting cavity mir- 
ror with a power transmission of 10 %. By operating the 
optical parametric oscillator below threshold, a squeezed 
vacuum state is generated^ One of the cavity mirrors 



is mounted on a piezo-electric element in order to adjust 
the optical path length of the cavity. The optical para- 
metric amplifier is stabilized by recording the transmit- 
tance through the cavity of a counter-propagating light 
beam on detector Dl (red dashed line) using the Pound- 
Drever-Hall technique that stabilizes the cavity over sev- 
eral hours^ A bright amplitude-squeezed state is gen- 
erated by overlapping the squeezed vacuum state with 
a coherent state on a non-polarizing beam splitter that 
transmits 99% of the squeezed vacuum state and mixes 
in 1 % of the strong coherent state, cf. Fig. [5h.. Note 
that the relative phase A<f> between the coherent state 
and vacuum squeezed state must be chosen properly in 
order to generate a bright amplitude-squeezed state. By 
continuously tuning the relative phase, Atf>, a light source 
with excess or reduced photon fluctuations relative to the 
classical limit is obtained. Both the pump beam and the 
coherent displacement beam are shaped in high finesse 
modecleaning cavities (MC) in order to ensure a good 
spatial mode matching of the vacuum squeezed state and 
the coherent state on the beam splitter. 

In order to characterize the bright squeezed light 
source, we use a highly sensitive InGaAs-resonance de- 
tector that converts the photo current of the photo diode 
(ETX 500T, Epitaxx) with a transimpedance amplifier 
into a voltage. An electrical spectrum analyzer computes 
the Fourier spectrum of the auto correlation function, 
T(r) (cf. Eq. flTJ), from the AC voltage of the photo 
detector at an electrical sideband frequency uj s with a 
resolution bandwidth of Suj/2tt = 300 kHz. The reso- 
nance frequency of our photo detector is 3.93 MHz, which 
is therefore the chosen sideband frequency u s /2-k — 
3.93 MHz in our experiment. Next, we relate the DC 
voltage of the photo detector to the mean photon flux 
by calibrating it with a power meter. To this end we 
measure the light power, P that is proportional to the 
mean photon flux and divide it by the energy of a sin- 
gle photon, i.e. {n(t)) = P/huj . The linearity of the 
photo detector is ensured by measuring the DC and AC 
voltage of the photodiode versus power of the coherent 
laser beam. The experiments are performed with optical 
powers of the light source larger than P = 5 /iW where 
the signal is 6.5 dB larger than the dark noise level of the 
detector. 

Figure^) shows the measured Fano factor, F(uj s ,Sui), 
of the light source as a function of the relative phase be- 
tween the coherent beam and the squeezed vacuum beam, 
A</>, as obtained by removing the multiple scattering sam- 
ple. The classical limit (gray triangles) is recorded by 
blocking the squeezed vacuum beam, thereby probing the 
bright coherent beam. The relative phase, A</>, is scanned 
with a mirror mounted on a piezo-electric element. We 
find that the photon statistics of the nonclassical light 
source can be continuously tuned below the classical limit 
(F(oj s ,5oj) = 1) to a minimum of F(u s , Slj) =0.52 ±0.02 
and above the classical limit to a maximum Fano factor 
of F(u s ,5lj) =4.6 ±0.2. 

We furthermore characterize the vacuum squeezed 
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state by the technique of quantum tomography where 
the Wigner function is reconstructed based on phase- 
sensitive measurements of the quantum fluctuations, see 
inset in Fig. Two photo detectors mounted at the 
output ports of a 50:50 beam splitter record the photo- 
current as a function of the time. Both time-resolved 
signals are subtracted from each other, mixed with an 
electronic local oscillator at a frequency of 3.93 MHz, and 
low-pass filtered with a bandwidth of 150 kHz. The rel- 
ative phase, Ac/), of the coherent light beam is slowly 
varied so that the Wigner function is reconstructed from 
the recorded data for different values of the phase, us- 
ing an iterative maximum-likelihood reconstruction algo- 
rithm. 2 — The Wigner function is a quasi-probability dis- 
tribution that contains the complete information about 
the quantum state. It is expressed in terms of the real 
and imaginary part of the quantized electric field, the so- 
called x and p quadratures. From Fig. [5£ we observe that 
the vacuum squeezed state has strongly reduced fluctu- 
ations in the x quadrature at the expense of enhanced 
fluctuations in the p quadrature fulfilling the Heisenberg 
uncertainty relation^ 

The multiple scattering experiment is conducted by fo- 
cusing the nonclassical light onto either the front surface 
or the back surface of the sample, cf. Fig. [5^, to perform 
total transmission and total reflection measurements, re- 
spectively. The multiply scattered light is collected with 
a microscope objective (NA — 0.63) and recorded with 
the photo detector D2. In the reflection geometry it 
is desirable to remove direct reflections due to single- 
scattering events from the multiply scattered distribution 
of light. To achieve this, the multiple scattering surface 
of the sample is illuminated under a steep angle of 69°, 
whereby the direct reflection does not pass through the 
collection microscope objective. 
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FIG. 3. a, Scanning electron microscope image of the multiple 
scattering medium consisting of HO2 that has been grinded 
into strongly scattering particles with a typical size of 200 nm 
and deposited on a glass substrate, b, Scanning of a sample 
at different positions with the tip of a micrometer screw to 
obtain the average sample thickness. The scanned positions 
are shown as squares and their color represents the measured 
thickness. The average thickness for this sample is determined 
to be L = (6 ± 2)/i m. c, Measured inverse total transmission 
through the multiple scattering samples versus sample thick- 
ness obtained using an integrating sphere. The line represents 
a linear fit to the data whereby the transport mean free path 
and the effective refractive index are obtained. 



B. Sample characterization 

The multiple scattering samples are prepared by 
spreading suspensions of titanium dioxide particles on 
microscope cover glasses resulting in a typical sample size 
of 45 mm x 25 mm. A scanning electron microscope im- 
age of a sample is shown in Fig. [3^- After evaporation of 
the liquid, the thicknesses of the samples are measured by 
scanning each sample using the tip of a micrometer screw. 
The statistics in the thickness determination is obtained 
by probing 20 different positions on the samples (Fig.[3h)- 
The central area of each sample is left untouched to avoid 
damage using this contact technique. Note that all thick- 
ness measurements are carried out far from the sample 
border where surface tension of the suspension will lead 
to large variations in the sample thickness. The aver- 
age thicknesses of the samples can be seen in Fig. [3b, 
where the error bars display the measured variance in 
the thickness. In the optical experiments discussed in 
the next section only a small region in the center of the 
sample is investigated, and we therefore expect that the 



recorded size fluctuations provide an upper bound for the 
uncertainty relevant in the optical experiment. 

The ensemble averaged total intensity transmission 
and reflection coefficients T, R are the characteristic pa- 
rameters of the multiple scattering medium in the ex- 
periments carried out here. Independent measurements 
of T are carried out by using an integrating sphere with 
two entrance ports with a detector mounted on one port 
and the samples on the second port. We record the total 
transmitted power, P, through a sample and the light 
power, Pq, without the sample to extract the total sam- 
ple transmission, T, and ensemble average by measuring 
at six different positions on the sample. Variations in 
the transmission coefficients were observed depending on 
which of the two sample surfaces were illuminated, since 
they are surrounded by different dielectric media (glass 
and air, respectively). The inverse of the total transmis- 
sion is in the thin sample approximation given by 30 
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Here z e \ and z e 2 are the extrapolation lengths of the 
medium that are related to the effective refractive index-2i 
and T sur f accounts for the reflections at the two sample 
boundaries that depend on the refractive indices of the 
surrounding materials (either n g i ass = 1.45 or ri a , r = 1.0). 
Fitting the experimental data of Fig. [3fc with this theory, 
we estimate a transport mean free path of I = 0.9±0.3 /im 
and an effective refractive index of n e g = 2.0 ± 0.4 corre- 
sponding to extrapolation lengths of z e i = 4.3 ± 0.3 fim 
and z e 2 = 4.6 ± 0.3 /im, respectively. From that the av- 
erage number of scattering events, N = (L cS /£) 2 , taking 
place in the multiple scattering media can be estimated, 
where L e g = L + z e i + z e 2 is the effective sample thickness 
including interface effects^ For our samples we estimate 
N « 300 — 1000 scattering events for sample thicknesses 
varying between L = 6 /im and L — 20 /im demonstrating 
that we are deep in the multiple scattering regime. 
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FIG. 4. Overall collection efficiency for each sample thickness 
as measured in transmission (t)t, gray triangles) and reflection 
(tir black circles). The error bars on the sample thickness are 
included only on the reflection data. 

A crucial issue for any quantum optics experiments 
in multiple scattering media is to collect a sufficiently 
amount of photons in order to see nonclassical features. 
Consequently, the total collection efficiencies in both 
transmission and reflection (j]t,r) are essential parame- 
ters that can be measured independently from the quan- 
tum noise measurements by recording the ratio of the 
total transmitted (reflected) power collected with the 
squeezing detector relative to the total power transmit- 
ted through (reflected from) the sample. In the reflection 
setup, the multiple scattering medium is illuminated un- 
der a steep angle of 69° . The fraction of the incident light 
that is multiply scattered is R = T sur {(a)R s . The light 
that enters the sample, T sux f(69 ) « 72%, is calculated 
from Fresnel relations at the sample surface. The amount 
that is reflected inside the sample, R s = 1 — T s , is ob- 
tained from sample transmission measurements with an 
integrating sphere using T — T s T sur f(0°) and incorporat- 
ing the surface reflection for light under normal incidence 
to the multiple scattering medium, see Eq. (fT4")l . The 
detection efficiency for reflection measurements is then 



given by i]n — (Pdct/-Po) /R, where Pq is the power of 
the light source and Pdct denotes the ensemble-averaged 
power of the reflected multiply scattered light detected 
by the squeezing detector. The measured detection effi- 
ciencies for the transmission and reflection experiments 
are presented in Fig. 2J They are found to vary with 
the thickness of the sample since the diffuse light exit- 
ing the random medium is collected with different effi- 
ciencies depending on the extent of the diffusion process. 
We note that the collection efficiency obtained with the 
microscope objective is much higher than that of an in- 
tegrating sphere (rj ~ 1 — 10%), although the numerical 
aperture (NA) of the objective limits the obtainable col- 
lection efficiency. 



C. Experimental results 
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FIG. 5. Measured Fano factor Ft,r(ui s , 8oj) after multiple 
scattering of a squeezed state, as recorded in a, reflection 
geometry and b, transmission geometry for L = 16 fim (black 
circles). The classical limit is depicted with gray triangles, 
c, Ensemble-averaged transmitted (filled gray symbols) and 
reflected (empty black symbols) Fano factor depending on 
the incident Fano factor, F. The inset shows a zoom in of the 
data that have been taken for Fano factors smaller than F < 1 
highlighting the quantum regime of multiple scattering. The 
sample thicknesses are 6 fim (triangles) and 20 /im (circles), 
respectively. The theory expressed by Eq. ([9| is shown by 
the straight lines and the classical limited is given by the 
dashed line. The ensemble average has been performed for 
six different sample positions and the result has been verified 
18 times per sample position. 

The Fano factor of multiply scattered light, mea- 
sured in reflection geometry and depending on the rel- 
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ative phase A</>, is displayed in Fig. [5^ for a sample 
of L = 16 fim. Using the bright squeezed light source, 
we observe photon fluctuations of the multiply scattered 
photons that clearly differ from the classical limit. Either 
classical or nonclassical photon fluctuations are recorded 
depending on Ac/). For comparison, Fig. [5b shows the 
Fano factor of the total transmission. In both cases the 
measurement is performed at a single sample position, 
representing a single realization of disorder. The ob- 
served variation in the Fano factor of the transmitted 
light is smaller than in the reflection. This can be ex- 
plained by a higher reflection than transmission of our 
sample since the Fano factor is proportional to the sam- 
ple transmission and reflection, respectively (cf. Eq. ((9|)). 

Multiple scattering is a random process and Fig. [5h., 
b are examples of data recovered for a single realization 
of disorder that consequently are unpredictable since all 
physical observables in multiple scattering are statisti- 
cal quantities. In order to compare the experiment with 
theory, we perform an ensemble average over different 
realizations of disorder. First, we tune the relative an- 
gle Acf> over 18 full periods (0, 367r) at a single sample 
position. Data taken at A</> + 2ttN for different N cor- 
respond to the same incident quantum state of light and 
can be used to average in order to account for fluctua- 
tions in the light source. The statistical ensemble average 
is obtained by repeating the procedure for six different 
sample positions. Figure |5fc shows the transmitted and 
reflected, ensemble- averaged Fano factor depending on 
the incident quantum state of light for two different sam- 
ples thicknesses. The output Fano factors are found to 
depend linearly on the Fano factor of the light source, 
F(ui s ,5ui), in very good agreement with the quantum 
theory of multiple scattering, cf. Eq. ©, after including 
the independently measured collection efficiency (Fig. 0} . 
For a nonclassical light source (F{ui s , Slo) < 1), the mul- 
tiply scattered photons have reduced photon fluctuations 
below the classical limit even after performing the ensem- 
ble average demonstrating that nonclassical properties of 
light can be detected after the random process of multi- 
ple scattering (inset Fig. [5]:). This phenomenon that is 
of purely quantum origin could open new possibilities to 
increase the information capacity in communication ap- 
plications since the amount of transmitted information 
can be increased in a strongly scattering environment* 3 ^. 
The classical limit of the information capacity is given 
by the signal-to-noisc ratio of the detected light utilizing 
a light source with Poissonian photon statistics. Using 
nonclassical light as in the present work breaks this clas- 
sical limit and therefore further increases the information 
capacity. 34 

Finally, we study the spatial correlations induced by 

multiple scattering of squeezed light. C® b , (uj s , Slu) does 
not depend on the spatial separation of the modes b and 
b' in the diffusive regime, i.e. the spatial quantum cor- 
relation is of infinite range. Therefore the correlation 
function can be determined from the photon fluctuations 
of the total transmittance or reflectance, as predicted 




1 2 3 4 5 

F(oo,8a)) 



FIG. 6. (color online). Measured spatial quantum correlation 
function depending on the incident Fano factor for different 
sample thicknesses (blue triangles: 6/j.m, black circles: 20/im). 
The spatial quantum correlation function has been measured 
in transmission (filled symbols) and reflection (empty sym- 
bols), respectively. The solid gray line displays the theoretical 
prediction while the dashed gray line represents the quantum 

limit of C?,, (u) s , Sco) = 0. For visibility we only show selected 
data points for each measurement. 



in Eq. (|13l) . Thus, we record the ensemble-averaged 
Fano factor Ft.r(w s , Suj) and the ensemble-averaged light 
power of the multiple scattered light. Using the cal- 
ibrated photo detector and an incident light power of 
P = 120 and hu = 1.17 eV (A = 1064 nm) we 
determine the incident mean photon flux get (n(t)) = 
P/hoj Q = 6.46 x lO^s" 1 . 

Figure [6] shows that the spatial quantum correlation 
function in a multiple scattering medium can be con- 
trolled by the quantum state of the incident light source. 
The correlation function is obtained from quantum noise 
measurements of the transmitted and reflected light, as 
outlined in Sec. IIB. The data are shown for two differ- 
ent samples {L = 6 /im and L = 20 /im) as a function of 
the incident Fano factor, F(ui s ,5uj). We clearly observe 

a linear dependence of C bb , (uj s ,8uj) on the Fano factor 
of the light entering the medium. The incident Fano 
factor and photon flux were measured by removing the 
sample. In excellent agreement with theory (Eq. (|12p). 
we fi nd in transmission and reflection measurements that 
C? b , (uj s , 8<jj) is independent of the sample thickness. The 
strength of the spatial quantum correlation varies over a 
broad range between —2 x 10~ 10 and 1.8 x 10~ 9 . Note 
that in Ref. El the bandwidth of the measurement was 
not accounted for properly, which underestimated the 
magnitude of the spatial correlation function 3 ^. We high- 
light that our observation of non- vanishing spatial quan- 
tum correlations that are tunable above and below the 
classical limit is in stark contrast to the classical intensity 
correlation function that is always unity in the diffusive 
regime, — The experiment demonstrates that the nonclas- 
sical photon statistics can survive multiple scattering and 
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persists even after ensemble averaging. An interesting 
next step will be to study experimentally the quantum 
interference of two or more independent quantum states 
in a multiple scattering medium?^ 

IV. CONCLUSION 

In conclusion, we have presented a continuous-mode 
quantum theory of multiple scattering in order to prop- 
erly describe spatial quantum correlations that are in- 
duced by multiple scattering. Depending on the quantum 
state of the incident light source we predict and observe 
that the multiple scattered light can exhibit either classi- 
cal or nonclassical photon fluctuations. From total trans- 
mission and total reflection measurements we observed 
either classical or nonclassical spatial correlations when 
varying the quantum state of the squeezed light source. 
Controlling the quantum state of the incident light source 
provides an efficient way of tuning continuously the mul- 



tiple scattering spatial quantum correlations. The ex- 
perimental results were in excellent agreement with the 
theoretical analysis. We expect our results might inspire 
new experiments on, e.g., multiple scattering of entan- 
gled photons or quantum interference! 12 ! 18 ' 19 ! 37 
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